Momentum distribution of confined bosons: temperature dependence 



J. Tempere 1 , F. Brosens 1 , L. F. Lemmens 2 , J. T. Devreese 1,2 '* 
1 Universiteit Antwerpen (UIA), Universiteitsplein 1, B2610 Antwerpen, 

Belgium. 

2 Universiteit Antwerpen (RUCA), Groenenborgerlaan 171, B2020 Antwerpen, 

Belgium. 

(Accepted in Phys. Rev. A on July 1, 1998) 

The momentum distribution function of a parabolically confined gas of bosons with harmonic 
interparticle interactions is derived. In the Bose-Einstein condensation region, this momentum 
distribution substantially deviates from a Maxwell-Boltzmann distribution. It is argued that the 
determination of the temperature of the boson gas from the Bose-Einstein momentum distribution 
function is more appropriate than the currently used fitting to the high momentum tail of the 
00 ' Maxwell-Boltzmann distribution. 
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I. INTRODUCTION 



The extreme low temperature required for the observation of Bose-Einstein condensation alkali vapors is a major 
achievement of present day experimentation [1-3]. Unfortunately, the temperature itself is not directly observable. 
' It has to be estimated from other, more accessible properties of the system, such as the velocity distribution of the 
\Q , vapor [4] . Because of the uncertainty relation between the inverse temperature and the energy of the system for a 
finite number of particles, it is not clear a priori how many particles are required for a reasonable thermodynamically 
defined concept of temperature. In the present paper, we do not address this fundamental question, but rather assume 
that the number of bosons in the system is sufficiently large to have a meaningful definition of the temperature. 

Under this assumption, the experimental estimate of the temperature [4] is based on the hypothesis that the 
momentum distribution function is Maxwell-Boltzmann like, at least for sufficiently high momentum. This hypothesis 
deserves further investigation and quantitative comparison to the momentum distribution of a parabolically confined 
Bose-Einstein gas. To the best of our knowledge, the Bose-Einstein momentum distribution in a parabolic confinement 
potential has not been documented before. 

Using the methods for path integration developed in [5-7] , we have been able to calculate the momentum distribution 
function of a gas of parabolically confined bosons with quadratic harmonic interparticle interactions. This calculation 
is presented in Sec. II. Other investigations using the same methodology [8] indicated that for weakly interacting 
Q H . systems, such as 87 Rb and 31 Na, the interaction slightly influences the parameters involved in the model system, and 
therefore have only a moderate influence the momentum distribution. 
• i-h , In section III, the approximations are examined which lead to the Maxwellian momentum distribution function, 
and a quantitative comparison is made between the temperature obtained in the Bose-Einstein distribution and the 
estimate of the temperature, obtained by fitting a Maxwell-Boltzmann distribution to the tail of the measured velocity 
distribution. The conclusion from this analysis is that the estimate from the Maxwellian deviates substantially from 
the temperature of the Bose-Einstein distribution in the presence of a Bose-Einstein condensate. 

II. MOMENTUM DISTRIBUTION FUNCTION OF PARABOLICALLY TRAPPED BOSONS 

The anisotropic one-body potential energy V\ and the two-body potential energy V2 of the model system consisting 
of N bosons are given by 

N 2 N 

Vi = ~ E Kv^ + ft + ^i) and V 2 = T— |— £ (*i r <) 2 W 
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where rj = (xj,yj, Zj) is the position of the jth boson, m is the mass of the bosons and £l xy , £l z are the frequencies 
of the parabolic confinement in the xy-plane and along the z-axis, respectively. The model interparticle interaction 
is harmonic with frequency u> and can be either attractive (-) or repulsive (+). The renormalized frequencies in 

the xy-plane and along the z-axis will be denoted by w xy = \jQ xy T Nlj 2 and w z = \j£l 2 =F Nuj 2 , respectively. 

These frequencies are the eigenfrcquencies of the internal degrees of freedom of the boson gas, whereas the confining 
frequencies fl xy , Q z are the cigcnfrequencics for the center-of-mass motion. The momentum distribution n DO sc(p), 
giving the statistical mean number of bosons with a given momentum p in the canonical ensemble, can be written as 



n hosc (p) = ( y${p- Pi) > (2) 




qp '|, (3) 

where the angular brackets indicate quantum averages for boson statistics. These mean values can be cal- 
culated for a system at temperature T = l/(fcs/3) using the path integral propagator in momentum space 

W,..,pMpi,..,p5 v ;0): 

" JV 



jPf/dpj A(p 1 , ..,pat)/C(pi, -,pjv;/3|pi, -,Pn;0) 



( Up, ..p.v); = VJ 1 /jy 7 r • (4) 

Jl/dpj /C(pi,..,pjv;/?|pi,..,pjv;0) 

The calculation of the momentum distribution, analogous to the calculation of the spatial density distribution n(r) 
in [7], can be summarized as follows. First the propagator in momentum space is evaluated by a Fourier transform of 
the corresponding propagator in position space obtained in [7]. This results in 

/C(pi,..,pjv;/3|pi,--,Piv;0) = /C d (pi, ..,pat;/3|Ppi, ..,Ppat;0), (5) 

permutations P 

with 

£d(pi, ..,pjv;/?|Ppi, ..,Ppjv;0) 



w z sinh {(3hw z ) w 2 xy sinh 2 {(3hw xy ) 
^ n z smh((3hn z )ni y smh 2 ((3hn xy ) jl^W'PW'") 

Pi /tanh(/3MV 2 ) tanh (/3?iw z /2) 

I hmN \ Q z w z / , 

' <,XP S P x + P y Aanh (phSl xy /2) _ tanh (phw x y /2) ^ ' ' : " ' 

HmN \ Q xy w xy 



where P = Pj i s tri e total momentum, and K<i(p";(3\p';0) is the one-particle propagator in momentum space 

for a particle in an anisotropic harmonic potential. Next, the sum over permutations appearing in the expression (4) 
for the mean value e~ lqp j ^ is transformed into a cyclic sum (with the number of cycles of length I). The 

cyclic summation can be performed for the generating function 

oo / JV \ 

N=l \j=i I 

where Z(N) = (rijLi / ^Pj) ^-(Pi> --iPn', P\pi, ■■, Pn', 0) is the partition sum. The result of the cyclic summation is 

hmq 2 ( Q z w z 



AN \coth(l3hn xy /2) coth((3hw xy /2) 
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u l exp {-jp [g> z coth(^w z /2) + (gg + gg)w z coth(l/?ftw z /2)] } 
8 sinh {l[3Tiw z /2) sinh 2 {£(3hw xy /2) ' 



(8) 



The expectation value e * q p J ^ for A^ bosons can then be derived in closed form from the generating functional 

gi(u) by collecting all the terms with the same power in u: 




-»q-p. 



exp 



n 2 



Timql 

AN y C oth(PhQ z /2) coth(/3hw z /2) 



hmjql + gj) 



AN 



xy 



W 



■nj 



x ^Z(N-t)^ 



coth(l3hn xy /2) coth(/3hw xy /2) J 

Km \ Qz w z cofh(£(3Tiw z /2) 
4 +(<lx + £) w z coth(lf3hw z /2) 



Z{N) 



8 sinh {£(3hw z /2) sinh'' {£(3hw xy /2) 



(9) 



Finally, the momentum distribution (2) is found as a Fourier transform of the expectation value (9). The result is: 



"•boso(p) = 



/ 1 



3/2 N 

£ 



Z(N - 1) 



\^hm) ^ Z{N) sinh 2 (SM^j ginh (^mj^ 

(pI+pD 



exp 



Pi 



hmAe(i} xy ,w X y) hmAe(n z ,w z ) 



At (Q xy , w xy ) y/A e (fl z ,W z ) 



with 



At (n, w) = 



w 1 

ta,nh(/3hw£/2) + N 



9. 



tanh(/37iO/2) tanh(/3?iw/2) 



(10) 



(11) 



In figure 1 this momentum distribution is shown at different temperatures for an isotropically confined Bose gas 
with w = fl and consisting of 1000 bosons. The momentum distribution for w/fl ^ 1 is qualitatively similar to the 
distribution with w/il = 1 shown in figure 1. In the inset, the value of the momentum distribution in the origin 
is shown as a function of the temperature. As the temperature is lowered below the condensation temperature, the 
momentum distribution becomes more pronounced in the origin: the average number of bosons with momentum zero 
rises sharply. 



III. APPROXIMATE VELOCITY DISTRIBUTIONS 



The momentum distribution (10) calculated in the previous section for a gas of parabolically confined bosons can 
be related to the Maxwell distribution through a series of approximations. First, in the case of a non- interacting 
(cj — > 0), parabolically confined boson gas the momentum distribution (10) becomes 



(pI + pD 



Z{N - I) exp ■ 



ideal 
"bosc 



(P) 



1 



AirTim 



3/2 N 

£ 



hmtt xy coth(/3hQ xy £/2) 

Pi 

TimVL z coth(/3M7 z ^/2) 



z(N)n xy sinh ((3hn xy e) ^n z sinh {jmij) 



(12) 



Furthermore, according to Fcynman [9], the Maxwell-Boltzman statistics can be recovered by considering only the 
cycles of length I = 1 in the cyclic summations of (10). This means that permutations involving cycles of length 
larger than 1 play no role any more. Hence the momentum distribution function for distinguishable non-interacting 
particles in a parabolic potential is: 



"dist al (P) 



tanh 2 (/3S,O a;2/ /2) tanh(/3M7 z /2) 
(irhm) 3 n 2 x v n z 



1/2 



cxp < 



{pI+pD 



hm£l xy coth(/3hQ, X y/2) 

Pi 

hm,n z coth(/3hQ z /2) 



(13) 
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If moreover the parabolic confinement were neglected, which corresponds to %Q, xy ,%Q, z <C 1//3 in (13), one finally 
would be left with the Maxwell-Boltzmann momentum distribution for free particles: 

n maxwe n(p)=^j ex P |-/3 j. (14) 

In experiments on Bosc-Einstcin condensed atomic vapors (e.g. [4]), the temperature of a cloud of parabolically 
trapped bosons was estimated by fitting the Maxwell-Boltzmann distribution (14) to the tail of the experimentally 
measured momentum distribution n mea sured(p)- This can be done by first choosing a sufficiently large threshold 
momentum p c such that atoms with a momentum larger than p c do not belong to the Bose-Einstein condensate, and 
then minimizing J p>p (fi maxwc ii (p) — ^measured (p)) dp with respect to the temperature T appearing in the Maxwell- 
Boltzmann distribution. The temperature obtained by this minimization will be denoted as Tvib and the procedure 
will be referred to as the "Maxwell-Boltzmann fit procedure" . 

Since the Maxwell-Boltzmann momentum distribution is a rather crude approximation to the momentum distribu- 
tion (12) of a gas of non-interacting, parabolically trapped bosons, we rather propose to determine the temperature 
by fitting n cxpor i mon t (p) to the momentum distribution (12) derived in the present paper. This is done by minimiz- 
ing J all (n^g 1 (p) — n moasuro d (p)) dp with respect to the temperature T appearing in the momentum distribution 

n bose (p) f° r non-interacting, parabolically trapped bosons. We will denote this procedure as the "Bose-Einstein fit 
procedure" and the resulting temperature as Tbe. 

In order to compare the Maxwell-Boltzmann fit procedure with the Bose-Einstein fit procedure, we applied the 
former on a momentum distribution simulated with H^g 1 (p) (« ".measured (p) for weakly interacting bosons), at a 

given temperature T BE . That is, j p>p ^ (n maxW cii (p,T M b) - (p,T BE )) 2 dp was minimized with respect to T M b, 

keeping Tbe and p c fixed. 

The result Tmb of the minimization is compared to Tbe in Figure 2, where the ratio (Tbe — Tvib) /Tbe is shown as a 
function of Tbe at a fixed threshold momentum p c = h^Jhmw. In the Bose-Einstein condensed phase the discrep ancy 
between T B e and T MB increases when the number of bosons increases, whereas for T B e > T = Tiw {/N/( (3), the 
discrepancy decreases when the number of bosons is augmented, as shown in the inset of figure 2. For temperatures well 
above the condensation temperature To, the ratio (Tbe — Tvib)/Tbe is small, and the Maxwell-Boltzmann fit procedure 
gives a reliable temperature estimate in this high-temperature regime. The difference between T B e and T MB becomes 
appreciable below the condensation temperature Tq. The present analysis therefore shows that (Tbe — Tmb) /Tbe is 
not negligible in the Bose-Einstein condensed phase and that the Bose-Einstein fit procedure is more appropriate than 
the Maxwell-Boltzmann fit procedure for the determination of the temperature. 

IV. CONCLUSION 

In this paper, we calculated the momentum distribution function (10) for a Bose gas in an anisotropic parabolic 
trapping potential, including quadratic interparticle interactions. In the absence of interparticle interactions, this 
momentum distribution simplifies to the momentum distribution function (12) of an ideal boson gas in a harmonic 
confining potential. The Maxwell-Boltzmann distribution (14) follows from (12) by the further simplifications of 
neglecting the boson statistics and the parabolic confinement. The underlying assumption in these derivations is 
that the number of bosons is sufficiently large to maintain the temperature concept from the thcrmodynamical limit. 
Assuming furthermore thermodynamical equilibrium, two temperature estimates were compared. Since the temper- 
ature estimates used in experiments on Bose-Einstein condensed atomic vapors (e.g. [4]), are obtained by fitting the 
Maxwell-Boltzmann distribution (14) to the tail of the experimentally measured momentum distribution ?i mcasurc d(p), 
we derived the corresponding temperature Tvib by fitting the large-momentum tail of the Maxwell-Boltzmann distri- 
bution to the large-momentum tail of the momentum distribution of the Bose gas (12). This temperature estimate was 
compared to the temperature Tbe in the momentum distribution of the boson gas. We find that Tvib substantially 
deviates from the correct temperature Tbe in the Bose-Einstein condensed phase, and hence that the Bosc-Einstcin fit 
procedure is more appropriate than the Maxwell-Boltzmann fit procedure for the determination of the temperature. 

ACKNOWLEDGMENTS 

Part of this work is performed in the framework of the FWO projects No. 1.5.729.94, 1.5.545.98, G. 0287.95, 
G. 0071. 98, and WO.073.94N (Wetenschappelijke Onderzoeksgemeenschap over "Laagdimensionele systemen", Sci- 
entific Research Community of the FWO on "Low Dimensional Systems"), the "Interuniversitaire Attractiepolen - 



4 



Belgische Staat, Diensten van de Eerste Minister - Wetenschappelijke, Technische en Culturele aangclegcnhcdcn" , 
and in the framework of the BOF NOI 1997 projects of the Universiteit Antwerpen. Two of the autors (J.T. and 
F.B.) acknowledge the FWO for financial support. 



[1] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Science 269, 198 (1995). 
[2] K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. Van Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, Phys. Rev. Lett. 
75, 3969 (1995). 

[3] C. C. Bradley, C. A. Sackett, J. J. Tollett, and R. G. Hulet, Phys. Rev. Lett. 75, 1687 (1995). 

[4] J. R. Ensher, D. S. Jin, M. R. Matthews, C. E. Wiemann, and E. A. Cornell, Phys. Rev. Lett. 77, 4984 (1996). 

[5] L. F. Lcmmens, F. Brosens, and J. T. Devreese, Phys. Rev. E 53, 4467 (1996). 

[6] F. Brosens, J. T. Devreese, and L. F. Lemmens, Phys. Rev. E 55, 227 (1997). 

[7] F. Brosens, J. T. Devreese, and L. F. Lemmens, Phys. Rev. E 55, 6795 (1997). 

[8] J. Tempere, F. Brosens, J. T. Devreese, and L. F. Lemmens, submitted to Phys. Rev. 

[9] R. P. Feynman, Statistical Mechanics, Addison- Wesley, Reading, MA (1972). 

Figure Captions 

Fig.l: In Figure 1 the average number n(p)/n(p = 0) of particles with a momentum p, is shown at several 
temperatures for a Bose gas in an isotropic parabolic confinement, as function of the momentum. Units are chosen 
such that h = m = Q = 1. In this set of units, the momentum is expressed in p = V hmfl. I n the inset, the 
momentum distribution at p = is shown as a function of temperature, with T = HQ/ks y/N/((3). For this figure 
we have chosen w = O : the momentum distribution for w/Q, ^ 1 is qualitatively similar to the distribution for 
w/fl = l. 

Fig. 2: The relative difference between the temperature Tmb found by fitting a Maxwell-Boltzmann distribution 
and the actual temperature Tbe is shown in this figure, representing (Tbe — ^mbV^be versus Tbe in units T = 
(Hfl/ ks) \/ N / ((3) . This is shown for 500, 1000 and 2000 non-interacting bosons. In the inset, the dependence of 
(Tbe — Tmb) /Tbe on the number of particles is shown. 
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Tempere, Brosens, Lemmens, Devreese, fig. 1 




Tempere, Brosens, Lemmens, Devreese, fig. 2 




